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Abstract 

We investigate a model in which tiny neutrino masses are generated at the two-loop level by using scalar 
leptoquark and diquark multiplets. The diquark can be singly produced at the LHC, and it can decay 
into a pair of leptoquarks through the lepton number violating interaction. Subsequent decays of the two 
leptoquarks can provide a clear signature of the lepton number violation, namely two QCD jets and a pair 
of same-signed charged leptons without missing energy. We show that the signal process is not suppressed 
while neutrino masses are appropriately suppressed. 
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I. INTRODUCTION 



The Standard Model (SM) gauge symmetry of elementary particles based on the SU(3)c x 
SU(2)x x U(l)y has been tested very accurately. On the other hand, the existence of the neutrino 
masses has been established [1-6]. This is clear evidence of the new physics beyond the SM because 
neutrinos are massless in the SM. Since neutrinos are electrically neutral, they can be Majorana 
particles unlike the other SM fermions Q] • The reason why neutrino masses are very different from 
those of the other SM fermions might be the Majorana property of neutrinos. 

The most familiar utilization of the Majorana property to generate tiny neutrino masses is the 
so-called Type-I Seesaw mechanism in which SU(2)i-singlet right-handed neutrinos mediate in the 
tree diagram jjj. Because of the suppression by mass scales of new heavy particles, naturally light 
neutrinos can arise. Another typical prescription to obtain tiny Majorana neutrino masses is the 



so-called radiative seesaw mechanism [91-4131] . where neutrino masses are induced at the loop level. 
In these models, the suppression of neutrino masses can be achieved by the loop suppression factor 
and/or a combination of new coupling constants even if new particles are not very heavy. The 
masses of charged leptons involved in the chirality flipping loop provide further suppression of the 
neutrino masses in some of such models [9- 11]. 

Although the lepton number is conserved in the SM, the addition of the Majorana mass term 
of neutrinos breaks the lepton number conservation by two units. The measurement of the lepton 
number violating (L#V) processes such as the neutrinoless double beta decay Q, is extremely 
important because it gives evidence that neutrinos are Majorana particles. Such processes are 
naively expected to be very rare because neutrino masses are very small. This is true for the Type- 
I seesaw model with very heavy right-handed neutrinos because light Majorana neutrino masses are 
unique lepton number breaking parameters at the energy scale which is experimentally accessible. 
However, in radiative seesaw models, a trilinear coupling constant for light (e.g. TeV-scale) scalars 
can be more fundamental than light neutrino masses as the L^V parameter at the accessible energy 
scale. Then, L^V processes via the trilinear coupling constant can be significant at the TeV-scale 
even if the neutrino masses are suppressed enough. 

New particles related to the neutrino mass generation are usually produced via the electroweak 
interaction, and therefore the production cross sections are not so significant at the LHC. However, 
new particles in the loop of the radiative seesaw models can be charged under the SU(3)c* [I6l4l8l|. 



Such a colored particle can easily be produced at hadron colliders. In these models, decay patterns 
of new colored particles could be related to the form of the neutrino mass matrix constrained by 
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TABLE I: List of particle contents of the model, 
the neutrino oscillation data |l6l4l9[|. 



In this paper, we investigate a radiative seesaw model with a scalar leptoquark multiplet and 
a scalar diquark multiplet. Major ana masses of neutrinos are induced via the two- loop diagram 
where colored particles are involved in the loop. The lepton number violation is caused by the 
trilinear coupling constant of the leptoquarks and diquark, which can produce a characteristic 
signature at the LHC. The signature consists of two QCD jets and a pair of same-signed charged 
leptons without missing energy, which would be easily observed at the LHC. 

This paper is organized as follows. In Sec. II, we present the model. Section III is devoted to 
discussion on the collider phenomenology and the low energy constraints for the leptoquark and 
the diquark in the model. Conclusions are given in Section IV. 



II. THE MODEL 



The particle contents of the colored radiative seesaw model are shown in Table HI The model 
is briefly mentioned in Ref. [if]]. The model includes a scalar leptoquark multiplet (Slq) whose 
lepton number and baryon number are 1 and 1/3, respectively. Under the SM gauge group, the 
Slq is assigned to the same representation of right-handed down-type quarks; a 3 representation 
of SU(3)c*, a singlet under SU(2)i, and hypercharge Y = —1/3. We also introduce a scalar diquark 
multiplet (Sdq) which has a baryon number 2/3. We take Sdq as a 6 representation of SU(3)c, a 
singlet under SU(2)^, and aF = —2/3 field. The diquark of a 6 representation can be expressed 
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in a symmetric matrix form as 

^ Sdqi 5 , D q4/v / 2 Sdqs/v^ 
Sdq= S mi /V2 S DQ2 S DQ6 /V2 ■ (1) 

y^DQs/V^ Sr>Qd/V^ SvQ 3 J 

The baryon number conservation is imposed to the model such that the proton decay is for- 
bidden. We introduce the soft-breaking term (see the next paragraph) of the lepton number 
conservation to the scalar potential in order to generate Majorana neutrino masses. The Yukawa 
interactions with the leptoquark and diquark, which preserve both of the lepton number and the 
baryon number, are given by 

^Yukawa = - {Zf < °2 Qf + Wf (Yr)u } (^Lq)* " J& (Xs)ij tf R )* + H.C., (2) 

where a a (a = 1-3) are the Pauli matrices, a and f3 (= r,g, b) denote the color indices; for example, 
5£q corresponds to 5dqi m Eq. (P). We choose the diagonal bases of mass matrices for the charged 
leptons and down-type quarks. Then, the SU(2)l partner of d iL is described as u' iL = (^ckm)jj u jL-> 
where Vckm is the Cabibbo-Kobayashi-Maskawa (CKM) matrix and Uj = (uj R ,Uj L ) T are mass 
eigenstates of up-type quarks. Mass eigenstates v iL of neutrinos are given by v iL = (ul ANS ) ii v iL , 
where C/mns is the Maki-Nakagawa-Sakata (MNS) matrix. The Yukawa matrices (Yj,, Yr, and Y s ) 
are 3x3 matrices under the lepton flavor {I = e, fi, r) and the down-type quark flavor (i, j = 1-3). 
While Yl and Y# are general complex matrices, Y s is a symmetric matrix (Yj = Y s ). Note that 
neutrinos interact with the leptoquark only through Yl, and we will see later that Yr is irrelevant 
to the neutrino mass at the leading order. 

In the scalar potential of the model, we introduce the following three-point interaction: 

f^(Sl Q T (S[ Q )*S^ Q + R.c. (3) 

The coupling constant ji softly breaks the lepton number conservation by two units while the baryon 
number is conserved. There is no other possible soft-breaking term of the lepton number and/or 
the baryon number. We can take the \i parameter as a real positive value by using the rephasing 
of S*dq- Considering radiative corrections to m r p and w, D q via the \x parameter, perturbativity 
requires fi < min(m L Q,m D Q) as discussed in Ref. 20 j for the Zee-Babu model (ZBM) [lo| . 

The neutrino mass term ^(M„)«/ VJ£ {i^^iY in the flavor basis is generated by a two-loop dia- 
gram in FIG. [T] including the leptoquark and the diquark. The mass matrix is calculated as 

(M„)«, = +24/X (Y* L ) H m d (Y s ) tj I tj m d (Y L %„ (4) 
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FIG. 1: The two-loop diagram for the neutrino mass generation in the model. 



where the loop function is defined as 
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The diagram is similar to the one in the ZBM although SU(3)c-singlet particles in the loop are re- 
placed with colored particles . Thus, we refer to this model as the colored Zee-Babu model (cZBM). 



(5) 



See, e.g., Refs. 



2014221] for studies about the ZBM for comparison with the cZBM. Note that Yr 



does not contribute to the two-loop diagram. 1 In the ZBM, at least one massless neutrino is pre- 
dicted because of the antisymmetric Yukawa coupling matrix. In contrast, all of three neutrino 
masses can be non-zero in the cZBM because Yl is not an antisymmetric matrix. Since new colored 
scalars should be much heavier than the SM fermions, the loop function can be reduced to 

1 1 



2l| 



Iij ~ J = 



(4tt) 4 (max[m LQ ,m DQ ]) 2 3 



(6) 



where 



l + -{(lnr) 



7T 



1} for r> 1 
for rCl 



(7) 



Hereafter, we restrict ourselves to the simplest scenario where Yr is small enough to be ignored. 2 
A benchmark point in the parameter space of the cZBM is shown in Appendix lAl 



1 The Yr contributes to Majorana neutrino masses at the higher loop level. The four-loop contribution is utilized 
in the model of Ref. where Yl is ignored. 

2 If we extend the model as a two-Higgs-doublet model, we can eliminate the Yr term by using a softly-broken Z2 
symmetry (e.g., u iR (or £r) and the second Higgs doublet are ^2-odd fields) which is also required to avoid the 
flavor changing neutral current at the tree level. Another example to eliminate the Yr term is the case where the 
leptoquark is not an SU(2)£-singlet but a triplet. 
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III. NEW COLORED SCALARS AT THE LHC 



A. Leptoquark 



The main production channel of leptoquarks at hadron colliders would be the pair-creation from 
gg and qq annihilation 



could also be possible 



at the leading order 



24j. The associated production of Slq with a lepton from qg coannihilation 



251 ] . The pair-production cross section is determined only by QCD interaction 



241 ] . while the associated production mechanism highly depends on the Yukawa 



coupling constant of the leptoquark 



251 ] . The leptoquarks have been searched at the Tevatron 



and the LHC. The most stringent lower bound on the leptoquark mass at 95% confidence level 
is set as 830 GeV (840 GeV) by the recent CMS result at *Js = 7TeV with 5.0 fb" 1 integrated 



luminosity 



2o ]; the pair-production of scalar leptoquarks is assumed as well as a hundred percent 



decay branching ratio into the first (second) generation quarks and leptons. See also Refs 
28 ] for the ATLAS results with 1.03 fb -1 integrated luminosity. The analysis of the decay into 
third generation fermions would be performed in near future. The search strategies for the third 
generation leptoquarks have been studied in Ref. [29]. 



The leptoquark induces various LFV processes. At the tree level, four-fermion operators 
(two left-handed leptons and two left-handed quarks) are generated by integrating leptoquarks 
out. The constraints on such operators have been extensively studied in Ref. [30f|. Tables 3, 
4, 12, and 13 in Ref. [30|] are relevant to the cZBM. Especially, operators (6^7 /i // L )(lZ7^7 M 'U i ) 
and ^J^ iiv i l L)^ , Ln^ s L> are strongly constrained by the /i-e conversion search and the K me- 
son decay measurement, respectively. For the benchmark point given in Appendix we have 
|(Y L ) e i(Y L *) M i|/(4 v / 2G F m2 Q ) = 6.1 x 10~ n and \(Y L ) n (Y*) e/2 \/(AV2G F m 2 LQ ) < 10^ 7 which sat- 
isfy constraints shown in Ref. [30J], where Gf = 1.17 x 10~ 5 GeV -2 . 

At the loop level, effects of leptoquarks on charged lepton transitions, i.e., li — > ij'j, have 



also been studied 



311 ]. Since we assume that S^q has the Yukawa interaction only with the left- 



handed quarks (namely Yr = 0), the contribution from the top quark loop does not give a large 



enhancement of li — > Ijj. 3 Then, the branching ratio of fj, — > &f is calculated as 



BR(// -»• &y) 



3a 



EM 



2567rG|mt Q 



Yl4 



Cfl 



(8) 



where a EM = 1/137. For example, a benchmark point shown in Appendix lAl gives BR^ — > ej) 



3 It is known that the similar process b — > (induced by the uncolored charged Higgs boson) in the Type-II 
two-Higgs-doublet model is enhanced by the top quark loop [32| . 



6 



6.5 x 10 13 which satisfies the current upper bound (2.4 x 10 12 at 90 % confidence level) in the 



MEG experiment 



3J. 



Since we take Yr = 0, the sign of the leptoquark contribution to the leptonic g — 2 cannot be 
changed. It is worth to mention that the contribution of the leptoquark has an appropriate sign 
(the plus sign) 4 to compensate the difference between the measured value and the SM prediction 
for the muon g — 2. The preferred size of Yl is (YlY^)^ ~ 1 for jti L q ~ 1 TeV. In order to satisfy 
LFV constraints with this size of Yl, a simple ansatz is that Yl is a diagonal matrix. Note that we 
must take care about the constraint on (z^/y^ ' m l)(^lTm s l) ( see Table 12 in Ref. [scj) even if Yl 
is diagonal; the constraint on (Y^iO^)/^ is difficult to be satisfied with (Y s )i2 < 1 because Yl is 
related to Y s through the neutrino mass matrix. We could not find any viable example of such a 
parameter set although it might exist with more complicated structures of Yl and Y s . 



B. Diquark 



At the LHC, the diquark Sj)q in the cZBM would be singly produced by the annihilation of two 
down-type quarks. 5 The single production mechanism has an advantage to search for the relatively 
heavy diquark due to the s-channel resonance. The single production cross section is determined 
by (Ys)n, which is evaluated in Ref. 



351 ] as a function of the diquark mass with a fixed Yukawa 
coupling constant. The (Y s )n in the cZBM is less constrained by the neutrino oscillation data 
because its contribution to neutrino masses is suppressed by rn^/m^Q. If we assume (Y s )n = 0.1 
and m D Q = 4 TeV, the single production cross section a(dd — > S*dq) is about 5fb at the LHC with 
y/s = 14 TeV 35J]. Note that the CMS experiment at y/s = 7 TeV with 1 fb" 1 integrated luminosity 



excludes diquark masses between 1 TeV and 3.52 TeV at 95 % confidence level by assuming the 
diquark decay into two QCD jets for the Eq diquark which couples with an up-type quark and a 
down-type quark [36J. See also Refs. [37h3S|. 



The diquark induces flavor changing neutral current processes in the down-type quark sector. 
Especially, it gives tree-level contributions to K°-K°, B^-B® and B®-B® mixings, resulting in strong 



constraints on Y s . By using the notations in Ref. 



401 ] . the benchmark point in Eqs. (lAlh gives 



C 



K - -(y/) 11 (y s ) 22 /(2m 2 3Q ; 



0. Similarly we have C B = +1.2 x 10" 13 GeV -2 and C B 



These values satisfy the constraints obtained in Ref. [40J (see also Refs. 



41 



0. 



4 For a heavy scalar <j> which interacts with fi L and a light fermion ip L as JT[{iIjl) c 4>, its contribution to the muon 
g — 2 has the plus sign if the electric charge of cj> is greater than —2/3. See also Ref. [34| . 
The diquark can also be created in pair via the gluon-gluon annihilation. 
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The diquark in the cZBM decays into not only a pair of the down-type quarks but also a pair 
of leptoquarks. The fraction of fermionic and bosonic decay modes is calculated as 



r(s DQ i -> s[ Q S£ Q ) r 



(9) 



This formula is the same for the other diquarks because of the SU(3)c symmetry. We focus on 
the case where the ratio in Eq. Q is less than about unity such that the branching ratio for 
Sdq — > SlqSlq becomes 0(10)%. Subsequently, 50% of each leptoquark decays into an up- 
type quark (a down-type quark) and a charged lepton (a neutrino). Then, the model provides 
a characteristic signature in FIG. [51 whose final state consists of two QCD jets and the same- 
signed charged lepton pair without missing energy. The decay chain of the diquark can be fully 
reconstructed at the LHC. This signature can be a smoking gun for the lepton number violation 
because no lepton number is taken away by invisible particles. There is no SM background in 
principle because the SM conserves the lepton number. It would be also very rare that the SM 
process mimics the signal process because the leptons in the signal events are too energetic to be 
produced in the SM process. 

It should be emphasized that the event rate of the process is not necessarily suppressed though 
the full process picks up all new coupling constants relevant to the small neutrino masses (namely, 
Y s , fj,, and Y£). One reason for that is because the process in FIG. [2] does not have suppressions 
with the two-loop factor l/(167r 2 ) 2 and with down-type quark masses, which are used for tiny 
neutrino masses. The other reason is that on-shell productions of a diquark and leptoquarks are 

r i 2 

utilized as a(dd -)• Sdq)BR(Sdq ->■ SlqSlq) Yle,i BR(Slq ->■ ^l u il) '■> even if a partial decay 
width is controlled by a small coupling constant (e.g., Slq ^L u iL via (Xl)u), its branching ratio 
becomes sizable when the total decay width is also controlled by small coupling constants. In this 
scenario, the cZBM seems the new physics model which is the most easily probed at the LHC and 
takes us to the top of the energy frontier. 

For the benchmark point shown in the Appendix [Aj the ratio in Eq. ([9]) is 0.18 for which 85 % 
of Sdq decays into SlqSlq. Then, 15% (7%) of Slq decays into a charm quark (a top quark) 
associated with an electron or a muon. Decays into an up quark are negligible for the bench mark 
point. The decay into a tau lepton might not be reliable because it gives missing neutrinos. Even 
if Slq decays into a top quark, hadronic decays (68 %) of from the top quark decay have no 
missing energy. As a result, the cross section for L^V events without missing energy is about 
0.18 fb at the LHC with \fs = 14TeV for the benchmark point. 
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FIG. 2: The same-signed charged lepton signature without missing energy at the LHC. 



In the energy scale which is much below the new scalar masses, the diagram in FIG. [2] becomes 
a dimension-9 operator of six fermions. Such L^V operators up to dimension-11 have been studied 
in Refs. [lfjjQ]- The dimension-9 operator in the cZBM is highly suppressed by the inverse power 
of mass scales of new colored particles as well as by new Yukawa coupling constants. Therefore the 
collider signature in FIG. [2] does not conflict with the stringent constraints from the other lepton 



number violating observable such as the neutrinoless double beta decay 



violating rare decays r ± -»• PM ± M ± (M = w,K) [44|, M ± -»• M'^i ± l' ± (M = B,K,D) [45[ and 



14] . and lepton number 



t -»• b£ + e' + W 



46] 



IV. CONCLUSIONS 



We have studied a model for neutrino mass generation with the scalar leptoquark Slq and 
the scalar diquark <S*dq- Tiny Majorana neutrino masses are induced at the two-loop level where 
the colored particles are involved in the loop. The trilinear scalar coupling constant between two 
leptoquarks and a diquark is the only parameter of the lepton number violation in this model. The 
diquark can be singly produced at the LHC via the resonance mechanism, and it can decay into 
a pair of leptoquarks through the lepton number violating coupling. The leptoquarks can further 
decay into a charged lepton and an up-type quark. Thus, the model gives a distinctive signature at 
the LHC, namely pp — > Sdq — > <S'lq<S'lq — > £~£'~jj without missing energy, which would be a clear 
evidence of the lepton number violation. We have shown that the lepton number violating process 
is not suppressed because of on-shell productions and decays of the diquark and the leptoquarks, 
while Majorana neutrino masses are highly suppressed. 
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Appendix A: A benchmark point 



Here, we show a benchmark point of the model: 

/ 8.1 x 10" 5 4.0 x 10" 2 -7.0 x 10~ 3 \ 



Y L 



-4.9 x 10" 5 5.3 x 10" 2 
3.1 x 10~ 5 -2.3 x 10" 2 

/ 1.0 x 10" 1 



4.4 x 10~ 2 
8.9 x 10~ 2 

\ 



(Ala) 



V 







H = 1 TeV, 



-1.2 x 10~ 3 

m LQ = 1 TeV ' m DQ 





1.2 x 10~ 3 
3.9 x 10~ 5 

4 TeV. 



(Alb) 
(Ale) 



We define an overall constant of neutrino masses as C = 24/i/o — 5.9 x 10 -6 GeV -1 . 

The neutrino mass matrix is diagonalized as ^mns^^mns = diag(mi, ni2e ia ^ , m3e* a 3i) with 
C/mns which can be parametrized as 



MNS 



(l o o\ 

C 23 S 2 3 
\0 -S23 C 23 / 



C13 si 3 e" 

10 

15 



is\ 



\-s 13 e ld c 13 J 



I 



C12 Su 
"Sl2 C12 

1 



(A2) 



where Cij (sij) denotes cosOij (sinfljj). We use the following values: sin 

2^23 = 1, sin 2 20i3 = 0.1, 

sin 2 26*12 = 0.87, 5 = 0, Ami, = 2.4 x 10" 3 eV 2 > 0, and Am|, = 7.6 x 10~ 5 eV 2 . Matrices Y L and 
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Y s in Eqs. (|A1|) are constructed by assuming the following structures: 



Y L = U ] 



(l 





MNS 



m 3 



m 2 + m 3 



m 2 



rri2 + m 3 



V 



(Ys)u 










m 2 



m 2 




m 3 



m 3 



m 2 + m 3 / 

\ 





V m 2m 3 



yzm s mbC 
m 3 — m 2 



0^ 
y 
^0 Zy 



(A3a) 



(A3b) 



It is easy to 
We use x = 
{Y s )u = 0.1 



yzm s mbC z 2 m 2 C J 
see that M v with these matrices results in m\ = x 2 m\C{Y s )n : a 2 i = 0j an d 031 = tt- 
10- 4 , y = 0.07, z = 0.1, m d = 5 x 10~ 3 GeV, m s = 0.1 GeV, m fo = 4.2 GeV and 
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